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The production of four jets in electron-positron annihilation allows for measuring the strong coupling 
and the underlying group structure of the strong interaction simultaneously. This requires next- 
to-leading order perturbative prediction for four-jet observables. In this paper we describe the 
theoretical formalism of such a calculation with sufficient details. We use the dipole method to 
construct a Monte Carlo program that can be used for calculating any four-jet observable at the 
next-to-leading order accuracy. As new results, we present the next-to-leading order prediction for 
the thrust minor, 2/4 and C parameter (at C > 0.75) four-jet shape variables and the four-jet rates 
with the Cambridge jet clustering algorithm. 



I. INTRODUCTION 

Electron-positron annihilation into hadrons is the 
cleanest process to test Quantum Chromodynamics 
(QCD) in high energy elementary particle reactions. 
In this process the initial state is completely known and 
there is a lot of quantities, for instance the total cross 
section and jet related correlations, that depend on the 
long distance properties of the theory very little. These 
quantities can be calculated in perturbative QCD as a 
function of a single parameter, the strong coupling. For 
this reason the various QCD tests at electron-positron 
colliders can be regarded as experiments for deter- 
mining ag. 

The other ingredient of QCD, that is in principle free, 
is the underlying gauge group. Although by now nobody 
questions that QCD is based upon SU(3) gauge theory, 
the "full" measurement of QCD, that is the simultaneous 
measurement of the strong coupling and the eigenvalues 
of the quadratic Casimirs of the underlying gauge theory, 
the Cf and Ca color charges, is not a purely academic 
exercise. The possible existence of light gluinos ||] in- 
fluences both the value of and the measured value of 
the color charges (or, assuming SU(3)c, the value of the 
light fermionic degrees of freedom Nf). Thus the only 
consistent framework to check whether the data favor or 
exclude the additional degrees of freedom is a simultane- 
ous fit of these parameters to data. 

In principle any observable depends on these basic pa- 
rameters. The sensitivity of a given observable on the 
color charges however, is influenced by the fact that in 
perturbation theory the three gluon coupling appears at 
tree level first for four-jet final states. In the total cross 
section and for three-jet like quantities the adjoint color 
charge appears only in the radiative corrections. There- 



fore, four-jet observables seem to be the best candidates 
to measure the color factors. Indeed, during the first 
phase of operation of the Large Electron Positron Col- 
lider (LEP) four-jet events were primarily used for mea- 
suring Cp and Ca These measurements however, were 
not complete in the sense mentioned above. The lack of 
knowledge about the perturbative prediction for four-jet 
observables at 0{al) prevented the experimental collab- 
orations from fixing the absolute normalization of the 
perturbative prediction, therefore, as could not be mea- 
sured using the same observables. 

Recently the next-to-leading order corrections to vari- 
ous four-jet observables have been calculated |1Q-M. In 
this article we give sufficient details of our work 
and present several new results for next-to-leading order 
predictions of four-jet observables that were not pub- 
lished before. The important development that made 
possible these calculations was that the one-loop ampli- 
tudes for the relevant QCD subprocesses, i.e. for e+e^ 
4 partons, became available. In Refs. (ijjl^ Campbell, 
Glover and Miller introduced FORTRAN programs that 
calculate the next-to-leading order squared matrix ele- 
ments of the e'^e~ ^ 7* — s- qqQQ and qqgg processes. 
In Refs. 1 19 2^ Bern, Dixon, Kosower and Wienzierl gave 



analytic formulas for the helicity amplitudes of the same 
processes with the e^e~ —> ^ A partons channel in- 
cluded as well. For the sake of completeness, in our work 
we use the amplitudes of Refs. for the loop correc- 
tions. Although the tree-level helicity amplitudes for the 
e+e^ —> 5 partons subprocesses had been known pl|], we 
calculated them anew and present the results in this arti- 
cle in terms of Weyl spinors conforming with the notation 
used for describing the one-loop helicity amplitudes . 
We also present the previously unpublished color linked 
helicity dependent Born matrix elements for the e'^e~ — > 
4 partons processes, which are needed also for the next- 
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to-leading order calculation of the three-jet production in 
deep inelastic scattering and for vector boson plus two-jet 
production in hadron collisions. 

In Sec. II we give details of the analytical and numer- 
ical calculation and describe how we parametrise our re- 
sults. In Sec. Ill we present the complete 0{al) predic- 
tions for the four-jet rates using the Durham algorithm 
p2[ and the Cambridge algorithm proposed recently ||2|] 
and make a comparison of these algorithms from the per- 
turbation theory point of view. We show the next-to- 
leading order results for the shape variables thrust minor 
and 2/4 often used in experimental analyses and for the C 
parameter distribution at C > 0.75. Sec. IV contains our 
conclusions. We present analytic results for the four- and 
five-parton tree-level helicity amplitudes in Appendix A, 
and perform the color summation in Appendix B. 



II. DETAILS OF THE CALCULATION 

A. Cancellation of infrared divergences 

It is well-known that the next-to-leading order correc- 
tion is a sum of two integrals, 



^NLO _ / j^NLO 



da^^^" = / da'^ + / da 



(1) 



where dcr^ is an exclusive cross section of five partons in 
the final state: 



(2) 



and d(T^ is the one-loop correction to the process with 
four partons in the final state: 



j^da^ = y drw < |xi-'°°p|2 > j4 



(3) 



— the real and virtual corrections. Although we specify 
our formulas to the case of four-jet calculation, one can 
use the formulas of this section to obtain the m-jet cross 
section by simply changing 4 (5) to m (m -I- 1). 

The two integrals on the right-hand side of Eq. (|l|) are 
separately divergent in d = 4 dimensions, but their sum 
is finite provided the jet function J„ defines an infrared 
safe quantity, which formally means that 



(4) 



in any case, where the five-parton and the four-parton 
configurations are kinematically degenerate. The pres- 
ence of singularities means that the separate pieces have 
to be regularized, and the divergences have to be can- 
celled. We use dimensional regularization in c? = 4 — 2e 
dimensions |24|, in which case the divergences are re- 
placed by double and single poles in e. We assume that 



ultraviolet renormalization of all Green functions to one- 
loop order has been carried out, so the poles are of in- 
frared origin. In order to obtain the finite sum, we use a 
slightly modified version of the dipole method of Catani 
and Seymour ||2^ that exposes the cancellation of the 
infrared singularities directly at the integrand level. 

The reason for modifying the original dipole formal- 
ism is numerical. The essence of the dipole method is to 
define a single subtraction term dcr^ — the dipole sub- 
traction — , that reguralizes the real correction in all of 
its singular (soft and coUinear) limits. Thus, the two sin- 
gular integrals in Eq. (|I]) are substituted by two finite 
ones: 



^NLO 



da- 



NLO 



where 



and 





^ [da« 

5 


-da^] 






M 



(5) 



(6) 



(7) 



There are many ways to define the subtraction term, but 
all must lead to the same finite next-to-leading order cor- 
rection. Since the virtual correction is not positive defi- 
nite, then depending on the size of the subtraction term it 
may happen that either dcr^'^^ or da^^'~' is not positive 
definite. From numerical point of view the best situa- 
tion is when both are positive definite, so that numerical 
cancellation of terms with opposite sign does not occur. 
We define the subtraction term as a function of a pa- 
rameter a e (0,1] which essentially controls the region of 
the five-parton phase space over which the subtraction is 
non-zero such that a = 1 means the full dipole subtrac- 
tion (see subsection II. B). By tuning the value of a, we 
can achieve that we add two positive definite integrals 
for almost all values of the observable to obtain the full 
correction. We use an a ~ 0.1, which is advantageous 
also for saving CPU time: The large number of dipole 
terms and their somewhat complicated analytic structure 
makes the evaluation of the subtraction term rather time 
consuming. Constraining the phase space over which the 
subtraction is zero we can speed up the program. 

In spite that the five-parton integral in Eq. (|^) is fi- 
nite, we introduce a very small cutoff in the phase space 
around the singular regions. Such a cutoff does not alter 
the value of the integral, but helps avoiding the cancel- 
lation of very large numbers that could lead to arbitrary 
values close to the singularity due to the finite machine 
precision. This cutoff is useful, but is also dangerous: 
if the subtraction is not correct, the five-parton integral 
becomes finite, but incorrect. The third advantage of us- 
ing the parameter a is that such errors can be spotted 
by varying the value of a and checking whether the full 
correction is independent of this parameter. 
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B. Dipole formulas for final state singularities 

In this subsection we recall those dipole factorization 
formulas that are relevant to our calculation. We do this 
only to the extent that we can define the simple modi- 
fication to the original formalism and the explicit cross 
section formulas of our calculation unambiguously. For 
further details we refer to the original work of Catani and 
Seymour |2|. 

In the dipole method the subtraction term is a sum of 
several dipole terms, 



^daA= ^ y'dr(5)i?,,,fej4, 



(8) 



where the dipole I'ij,^ is a function of the final state 
momenta pi and is given by 



25y,fc(pi, ...,P5)=- 



1 



2pi ■ pj 



(9) 



X4(. ..,ij,...,k,...\ '^'^ Vy.fc |. . . , ij, . . . , fc, . . .)4 . 

In Eq. \. . . . . . ,k, . . .)i is a vector in color + he- 

licity space defining the four-parton amplitude obtained 
from the original five-parton Born amplitude by replac- 
ing a) the partons i and j with a single parton ij {the 
emitter) and h) the parton k with the parton k {the spec- 
tator). The momenta of the spectator and the emitter 
are given in terms of a dimensionless variable 



Pipj 



Pk 



1 



yzj,k 



Pk 



pipj + Pjpk + Pkpi 



Pi, ^Pt +Pi 



1 



Pk 

Uij.k 



(10) 



(11) 



Tk and are the color charges of the spectator and the 
emitter. These color charges are defined by their action 
onto the color space: If particle i emits a gluon with 
color index c then the color-charge operator has the 
following matrix element in color space 



[ci, ■ 



where T': 



, , . . . , C4, c\Ti\bi, .. .,bi,. 

hi ■ ■ ■ "0464 J 



,b4) = 



(12) 



Scibi ■ ■ ■ TZh. . ■ .5, 



: F^ f^, = —ifccibi (color-charge matrix in the 
adjoint representation) if the emitting particle z is a gluon 
and T^.h- = ^cib; (color-charge matrix in the fundamental 
representation) if the emitting particle z is a quark (in 



the case of an antiquark emitter T^.fj. = t 



In Eq. (|9|) the splitting matrices, Vij^k are matrices in 
the helicity space of the emitter. They depend on the 
kinematic variables ytj^k and 



PiPk 
PijPk 



(13) 



and take different forms for the splitting of different par- 



tons (see Eqs. (5.7-5.9) in Ref. l25|). 

The definition of the dipole momenta makes possible 
the exact factorization of the five-particle phase space 
into a four-particle and a one-particle phase space 



dr(5)(pi,.. 
drW(. 



where 



,P5) = 

• J Pij J 

d^p, 



(14) 



,Pk, ■ ■ .) [dpi{p^j,pk)] 



[<iPiiPi3,Pk)] j^^p^S+{p,)J{pf,p,j,pk) , (15) 
and the Jacobian factor is 
J{p^;P^J,Pk) = 6(1 - z,) 6(1 - y,,,k) ^^"y^'' ' . 

(16) 

As mentioned before, we modify the original formalism 
such that we constrain the phase space over which the 
dipoles are subtracted: 

Jda^{a) = /dr^'^ ^y.fc ©(y^.fc < a) ^4 , (17) 

with a € (0,1]. The jet function in Eq. ( pT[ ) is a four- 
parton jet function that depends on the momenta pij 
and Pk, but not on pi, therefore, the integral over the 
one-parton phase space can be performed analytically. 
It can be shown that after integration of the dipole 
2?ii,fc(pi, ■ • ■ ,P5) over [dpi{pij ,pk)], only color correla- 
tions survive |25[| , in the form: 

J [dp^iPij,Pk)] 'Dij^kipi, ■ ■ ■ ,P5) 6(yij,fe < a) = (18) 



Vy,fc(a) 



Mt'{ 



, k, . 



where 



■A) 



= 41, 



,4|T,;.T,|1,...,4)4 (19) 



are the color-correlated four-parton tree matrix elements. 
Their explicit expressions are given in Appendix B. In 
Eq. (H) 



V»j,fc(a) 



(20) 



[dpi(pij,Pfc)] Q{yij,k < a) 7^ < Vij^k > = 

^Pi ' Pj 

as 1 ( 4V ^ 
27rr(l-e) W,Pk) 

where < Vij^k > denotes the average of Vij,/c over 
the polarisations of the emitter parton ij. The func- 
tions Vij{e,a) depend only on the flavour indices i and 
j. Rewriting the one-particle phase space in terms of 
the kinematic variables Zi and y, from the definition of 
Vij (e, a) in Eq. (|20|) we obtain 
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Vy(e,a)= / dzi {zi{l - Zi)y 



(21) 



where the spin-averaged spUtting functions are given in Eqs. (5.29-5.31) of Rcf. |25j. Performing the integration in 
Eq. (|l|), we find 





a) 






a) 






a) 


= 2Ca 



1 ,2 



- ( 1 + a — log a 



16 



11 / I 



Vggie, a) = 2CAU^- log a ] + — [ - ~ I + a -loga ] + — ~ — + 0(e) 



6 



5-- + 0(e) 



50 TT^ 



9 



(22) 
(23) 
(24) 



It was shown in Ref. pq] that after integrating over 
the factorized one-particle phase space, the subtraction 
term can be recast in the form 



da' 



drW4(l,...,4| J(£,a) |1,...,4)4 J4, (25) 



where the insertion operator /(e, a) depends on the color 
charges and momenta of the four final-state partons in 
|1,...,4)4: 



-f(Pi,P2,P3,P4;£,a) 



1 



27rr(l-e) 



(26) 



The singular factors Vi(e,a), are defined as 

Vg{e,a) = ^Vgg{e,a) +N[Vgg{e,a) . 



(27) 
(28) 



Using Eqs. (22-^4|) they can be written in the following 
explicit form: 



'\e^ 3 

+0(£) , 



log" a] +ji ( - + a — log a 



(29) 

+ K, 



where the 7.^ and Ki constants are defined by 

3 ^7 



79(9) = 2'^'^ 



^9(9) = 



Ig^^^A-^TRNf, Kg-. 



67 
18 



Cf 



Ca 



(30) 



10, 



-TrNi 



The formal result of the cancellation mechanism dis- 
cussed in this subsection is that the next-to-leading or- 
der correction is a sum of two finite integrals as given in 
Eq. (||). We would like to mention that although both 
integrals are finite, the integrand of the five-parton in- 
tegral is in fact divergent, it contains integrable square- 
root singularities in the kinematically degenerate region 



of the five-parton phase space. The efficient way to in- 
tegrate such a function is to apply important sampling. 
We apply multichannel Monte Carlo integration for this 
purpose, but do not consider the details of this technical 
point in this article. 



C. The general structure of the results 

Once the phase space integrations in Eq. (|^) are carried 
out, the next-to-leading order differential cross section 
for the four-jet observable O4 at a fixed scale Q takes the 
general form 



1 da 
ao dOi 



(31) 



where 77 = as{Q) Cf/ 27r. The renormalization scale de- 
pendence of the cross section is obtained by the sub- 
stitution 77 — > ?7(/i) (1 -I- (3q/Cf ^'O-Xfi), with 7y(/x) = 
as{^Ji) Cf/ 27r, which yields 



1 do- 
fjo d04 



{Oi)^il{tifBoAOi) 



(32) 



BoAOi)^\nxl + CoAOi) 
Cf 



In Eq. (32) ctq denotes the Born cross section for the 
process e+e^ — s- qq, fj, is the renormalization scale, 
Xfi — mA/s is the renormalization scale divided by the 
total cm. energy and B04 and C04 are scale indepen- 
dent functions, -B04 is the Born approximation and C04 
is the radiative correction. We use the two-loop expres- 
sion for the running coupling. 



with 



asjMz) 
w{p,Mz) 



1 



I3i as{Mz)Hw{p,Mz)) 



2tt 



w{n,Mz) 



as(go) 
27r 



ln(^ 



(33) 



(34) 



4 
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B - 

Pi - 



(35) 



(36) 



with the normalization Tr. = 1/2 in Tr(T°Tt'') = ^R(5'^^ 
The numerical values presented in this letter were ob- 
tained at the Z° peak with Mz = 91.187GeV, Tz = 
2.49GeV, sin^ = 0.23, as{Mz) = 0.118 and iVf = 5 
light quark flavors. 

In order to make possible the measurement of the color 
factors, we write both the Born approximation and the 
higher order correction as linear and quadratic forms of 
ratios of the color charges [Ea : 



and 



where 



B4 ^ Bo + Bx X + By y , 

C4 = Co + C^x + Cyy + CzZ 
+ Cxx X + Cxy X y Cyy y 



(37) 
(38) 



(39) 



x = T^ , y = -Fr 

At next-to-leading order the ratio z appears that is re- 
lated to the square of a cubic Casimir, 

Na 

C3= ^ Tr(T"r''rt")Tr(Tt'=r''r") , (40) 

a,b,c—l 

C3 



via z 



The Born and correction functions Bj 



and Ci are independent of the underlying gauge group. 
In the next section we present the ,64 and C4 functions 
for various four-jet observables. 



III. RESULTS 

Four-jet observables can be classified into three major 
groups: (i) four-jet rates; (ii) four-jet event shapes; (iii) 
four-jet angular correlations. Detailed results for observ- 
ables falling into all three classes were already presented 
in the literature. Dixon and Signer gave full account of 
the next-to-leading order four-jet rates with three differ- 
ent (EO, Durham and Geneva) jet algorithms |ll[]. In 
Ref. [|l3| we confirmed their results. Among the four-jet 
event shapes the D parameter, acoplanarity, and the Fox- 
Wolfram moments Hi and 114 were calculated at 0{al) 
in refs. ||l^ an d [T^ ]. The results for the D parameter 
were confirmed pWT As for angular correlations Signer 
presented the leading color corrections in Ref. and 
we added the full corrections in Ref. . In this section 
we would like to add the four-jet rate obtained using the 
Cambridge clustering and several event shape variables 
to the list of four-jet observables that are calculated at 
the next-to-leading order accuracy. We do not consider 
the four-jet angular correlations here. 



A. Four-jet rates 

The most important multi-jet observables that are 
used for determining the underlying parton structure of 
hadronic events are the multi-jet rates. In e+e^ anni- 
hilation the widely known Durham algorithm have 
become indispensable for this purpose. Recently a new 
jet clustering, the Cambridge algorithm was proposed as 
an improved version of the Durham scheme [ p3[ . This 
scheme is designed to minimize the formation of "junk 
jets" — jets formed from hadrons of low transverse mo- 
menta, unconnected to the underlying parton structure. 
As a result, the hadronization corrections to the mean jet 
multiplicities were found smaller when the Cambridge 
algorithm is employed than for the Durham clustering 
p3| . However, it was shown in Ref. |^ that the small 
hadronization corrections found for the Cambridge algo- 
rithm in the study of the mean jet rate are due to can- 
cellations among corrections for the individual jet pro- 
duction rates. Apart from the very small values of the 
resolution parameter, j/cut < 10~'^'^, for the individual 
rates the Durham clustering shows comparably small (for 
2/cut > 10~^), or even much smaller hadronization correc- 
tions. In this subsection we present the next-to- leading 
order production rates for four jets using both algorithms 
and compare the size of the radiative corrections. 

The four-jet rates are defined as the ratio of the four-jet 
cross section to the total hadronic cross section: 



i?4 



g"4-jct 
Ctot 



(yc 



(41) 



V Bi{ycnt) + 7? ( C4(ycut) - -^B^iy cnt] 



where we used CTtot = cro(l + f??)- Setting the color 
charges to the SU(3) values, we plot the scale indepen- 
dent B4{ycut) and C4(2/cut) functions in Figs. 1 and 2 and 
tabulate the values for C4(?/cut) in Table I. 




-2.5 -2.0 

logio(ycut) 

FIG. 1. The Born function B4 for the four-jet rate as a 
function of the resolution variable j/cut with Durham (solid) 
and Cambridge (dashed) algorithms. 
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FIG. 2. The correction function C4 for the four-jet rate as 
a function of the resolution variable j/cut with Durham (solid) 
and Cambridge (dashed) algorithms. 



TABLE I. Correction functions to the four-jet rates for 
Durham and Cambridge algorithms. 



iogio(yc«t) 
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Comparing the values for the two Born functions, we 
see that at leading order the Cambridge algorithm gives 
slightly higher rates and the difference increases with de- 
creasing ?/cut- On the other hand, the correction functions 
become smaller for Cambridge clustering with decreasing 
j/cut- The result of these opposite trends is that the K 
factors, defined as 




-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 
logio(ycut) 

FIG. 3. K factors as a function of the resolution variable 
j/cut for Durham (solid) and Cambridge (dashed) algorithms. 

i^(2/cut) = l±r;(^/i)f^^, (42) 

are smaller for the Cambridge algorithm for small values 
of ycut, which is demonstrated in Fig. 3. 

The smaller K factors also mean smaller renormaliza- 
tion scheme dependence, which can be seen from com- 
paring Figs. 4 and 5. The usual interpretation of the 
smaller scale dependence is that the effect of the uncalcu- 
lated higher orders are expected to be smaller in the case 
of Cambridge clustering. It is interesting to note that 
in the middle j/cut region (lO"^-^ < j/cut < 10"^), where 
the hadronization corrections for the Cambridge cluster- 
ing were found significantly larger than for the Durham 
algorithm, the theoretical uncertainty due to the renor- 
malization scale ambiguity is smaller for the Cambridge 
than that for the Durham clustering. Of course, one 
has to keep in mind that the //-dependence bands are 
not upper bounds on errors that arise from truncation of 
the perturbation series, just suggestions. In particular, 
if there is an artificial narrowing of the /^-dependence 
bands, e.g. at a crossover point, they almost certainly 
do not represent the size of the truncation error at that 
point. 

Four-jet fractions decrease very rapidly with increasing 
resolution parameter ?/cut- As a result, most of the avail- 
able four-jet data are below ?/cut = 0.01. It is well-known 
that for small values of ^cut the fixed order perturbative 
prediction is not reliable, because the expansion param- 
eter as In^ Ucnt logarithmically enhances the higher order 
corrections. One has to perform the all order resum- 
mation of the leading and next-to-leading logarithmic 
(NLL) contributions. This resummation is possible for 
the Durham algorithm using the coherent branching for- 
malism ||2^ and the procedure is the same for the Cam- 
bridge algorithm |^^. The four-jet rate in the next-to- 
leading logarithmic approximation is given by p9| 
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FIG. 4. The QCD prediction for the four-jet rate with 
Durham clustering at Born level (light gray band) and at 
next-to-leading order (dark band). The two narrow bands 
show the four-jet rate in the NLL approximation {K — 0, 
lower band) and in improved NLL approximation (upper 
band) as explained in the text. The bands indicate the theo- 
retical uncertainty due to the variation of the renormalization 
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FIG. 5. The QCD prediction for the four-jet rate with 
Cambridge clustering at Born level (light gray band) and at 
next-to-leading order (dark band). The two narrow bands 
show the four-jet rate in the NLL approximation {K = 0, 
lower band) and in improved NLL approximation (upper 
band) as explained in the text. The bands indicate the theo- 
retical uncertainty due to the variation of the renormalization 



/%r,(Q,g)A<,(g,Qo) ) (43) 
JQo J 

dq rq{Q,q) Ag{q, Qq) 

dq' (Tgiq, q') A,(g', Qo) + (g') A/ (g', Qo)) 



In Eq. ( |43| ) the functions Aa{Q, Qo) are the Sudakov form 
factors which express the probability of parton branching 
evolution from scale Qo — Qy/Vcut to scale Q without 
resolvable branching. The Sudakov factors are defined in 
terms of the Pab{cts{q), z) vertex probabilities as follows 



Aa(Q,(3o; 
exp - 



(44) 



It was shown in Ref. |30|] that one can obtain an improved 
theoretical prediction for the differential two-jet rate if 
the vertex probabilities are taken at next-to-leading order 
tell, which we also consider in our analysis: 



l-z In 1 



Pqqi^s, Z) = Cp 

Pggias, z) ^ TRNf {z^ + (1 - zY' 



(45) 



^ — z ,^ , as 2 



The K coefficient is renormalization scheme depen- 
dent. In the MS scheme it is given by 



K^Ca 



67 
18 



TT 

~6 



10 



9 



TRNf 



(46) 



Performing the z integral in Eq. (^), one obtains the 
Sudakov factors as integrals of the emission probabilities 
Ta{Q, q) in the following form: 



A,(Q, Qo) = exp J^dqTgiQ, q)j 
A<,(Q,Qo) = 

exp ( - / dq [rg{Q,q)+Tf{q)] 



A/(Q,Qo) 



[\iQ,Qo)r 



As(Q,Qo) ' 
and the NLL emission probabilities are 



r,(Q,g) = 

TT q 
2Ca a,{q) 



"^('?).^Vn^ 



r/(0,9) 



1 

TT q 

Nfa,{q) 



2tt 
as{q) 



„ K In 
271 J q 



3 

4 

11 
12 



37r q 



(47) 
(48) 

(49) 

(50) 
(51) 

(52) 
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We relate the as{q) strong coupling appearing in the 
emission probabilities to the strong coupling at the rele- 
vant renormalization scale, according to the one- 
loop formula 



Durham algorithm 



(53) 



where w{q,qo) was defined in Eq. (34), and we use 
Eq. ( |33| ) for expressing as(M) i"^ terms of as{Mz) = 
0.118. We could also use a two- loop formula for as{q), 
but the result would differ only in subleading logarithms. 

The result of this resummation together with its renor- 
malization scale dependence is also shown in Figs. 4 and 
5 (narrow bands). The lower band corresponds to the 
usual NLL approximation (K = 0), and the upper band 
is the result of the improved resummation. We can see 
clearly from the figures that the fixed-order and the NLL 
approximations differ significantly. One expects that for 
large values of j/cut the former, and for small values of 
?/cut the latter is the reliable description, therefore, the 
two results have to matched. 

The Durham and Cambridge four-jet rates can be re- 
summed at leading and next-to-leading logarithmic or- 
der, but they do not satisfy a simple exponentiation 
psf . For observables that do not exponentiate the viable 
matching schemes are the R matching or the modified 
R matching |^,^ . We use R matching according to the 
following formula: 



r>R— match nNLL 



f (i34-Br^) 



(54) 



C. 



NLL 



(54 - ^4'''^'^) 



where B^^^ and C^^^ are the coefficients in the expan- 
sion of as in Eq. (|l|). 

In Fig. 6 we show the theoretical prediction at 
the various levels of approximation: in fixed order 
perturbation theory at Born level (LO), at next-to- 
leading order (NLO), resummed and R-matched pre- 
diction (NLO-I-NLL) and improved resummed and R- 
matched prediction (NLO-I-NLL+K). Also shown is the 
four-jet rate measured by the ALEPH collaboration at 
the peak corrected to parton level using the 
PYTHIA Monte Carlo H]. We used bin- by- bin correc- 
tion and the consistency of the correction was checked 
by using the HERWIG Monte Carlo |§. The two pro- 
grams gave the same correction factor within statistical 
error. The errors of the data are the scaled errors of 
the published hadron level data, and we did not include 
any systematic error due to the hadron to parton correc- 
tion. In the inset we indicated the renormalization scale 
dependence of the 'NLO-fNLL+K' prediction. 
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FIG. 6. The QCD prediction for the four-jet rate with 
Durham clustering in fixed order perturbation theory at lead- 
ing (dotted) and next-to-leading order (dashed), and fixed 
order matched with resummed (dashed-dotted) and improved 
resummed (solid) calculation compared to ALEPH data ob- 
tained at the Z'^ peak and corrected to parton level (error- 
bars). The renormalization scale is set to a;^ = 1. The lower 
plot shows the relative difference 84, — (data-theory) /theory, 
where theory means the next-to-leading order prediction 
matched with improved resummed calculation at a;,j = 1. The 
inset shows the renormalization scale dependence of the 'the- 
ory' prediction with scale variation 0.5 < < 2. 

Fig. 6 deserves several remarks. First of all, we see 
that the inclusion of the radiative corrections improves 
the fixed order description of the data using the nat- 
ural scale Xp, = 1 for larger values of ?/cut- Secondly, 
the importance of resummation in the small j/cut region 
is clearly seen, but it is still not sufficient to describe 
the data at the natural scale, neglected subleading terms 
are still important.^ On the other hand, the improved 
resummation seems to take into account just the right 
amount of subleading terms and it makes the agreement 
between data and theory almost perfect over the whole 
2/cut region as can be seen from the lower plot. (Although 
for 2/cut > 10~^'^ ^4 falls outside the ±5 % band, one 
should keep in mind that in this region the error of the 
hadron to parton correction is very large. Also, for the 
'NLO+NLL+K' prediction we found remarkably small 
scale dependence for ?/cut > lO^'^. This feature, however. 



* Our 'NLO+NLL' resuhs differ from those in Ref 0, 
where as{q) in calculating R^^^ was kept at the fix as{Mz) 
value EtII. 
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should be taken with care. The improvement, obtained 
by including the two-loop coefficient K, affects NNLL 
terms, but there are other contributions of the same or- 
der that are not taken into account (e.g., next-to- leading 
order running of as and other dynamical effects), which 
is not the case for the 2-jet rate. The scale dependence 
of the 'NLO-I-NLL-I-K' result would consistently be un- 
der control only after the inclusion of the complete set of 
NNLL terms. 

Finally it is worth noting that for j/cut — 10^^'^ both 
PYTHIA and HERWIG yield less than 2 % hadronization 
correction. At the same value of the resolution parame- 
ter the theoretical prediction is insensitive to corrections 
beyond next-to-leading order (the NLO, NLO-fNLL, 
NLO-I-NLL+K curves cross, the renormalization scale de- 
pendence is small), therefore, at this accidental value of 
Ucut the next-to-leading order prediction agrees perfectly 
with the hadron level data. 



B. Four-jet event shapes 

Four-jet event shapes were used extensively by the LEP 
collaborations for QCD studies 



37 1 . In this subsection 



we consider four shape variables, the 1/4 distributions for 
the Durham and Cambridge algorithms, the thrust mi- 
nor (Tinin) and the C parameter for C values above 0.75, 
which are often used in the experimental analyses. 

In the case of event shape distributions we multiply 
the normalized cross section with the value of the event 
shape parameter, so we use the paramctrisation 



E(04)^04-;^(04) 

(Jo dC4 



B{0,)^^\uxl 



(55) 
C{Oi) 



instead of Eq. (32), in which case the average value of 
the shape variable is easily obtained from the differential 
distribution: 



< O4 >5= d04S(04) 



(56) 



Using this paramctrisation we define the K factors of the 
differential distribution as 



KiOi) = 1 + i^i^fs) 



CjOj) 
B{Oi) 



(57) 



In the following wc plot the physical cross sections £(04), 
the K{04) factors and tabulate the correction functions 
C{Oi) for O4 = 2/4, T^in and C. 

The 2/4 value denotes the transition value for j/cut at 
which, when decreasing j/cut, the classification of a given 
event changes from three jets to four jets. The advantage 
of this variable over the differential four-jet rate is that 



this variable can be defined on an event by event basis. 
Depending on the actual resolution variable one obtains 
the distribution for the Durham clustering and the 
distribution for the Cambridge clustering. We calculated 
the B{yi) and C(j/4) functions for both algorithms. The 
B{yi) values equal the ycutB{ycut) values when y^ = ^cut, 
therefore, we tabulate only the C^y^) functions for the 
two algorithms in Table II. 

TABLE II. Correction functions to the differential distri- 
butions of the 1/4 variables for the Durham and Cambridge 
algorithm. The parameter values are at the lower edge of the 
corresponding histogram bin. 













0.000 


(2.523 ±0.425) 


• 10^ 


(1.064 ± 0.350) 


• 10^ 


0.005 


(2.212 ±0.017) 


• 10^ 


(1.857 ±0.019) 


• 10^ 


0.010 


(1.376 ±0.009) 


• 10^ 


(1.166 ±0.011) 


• 10^ 


0.015 


(9.429 ±0.071) 


■ 10^ 


(8.144 ±0.080) 


• 10^ 


0.020 


(6.799 ±0.070) 


• 10^ 


(5.855 ± 0.062) 


• 10^ 


0.025 


(4.930 ±0.063) 


• 10^ 


(4.346 ± 0.051) 


• 10^ 


0.030 


(3.760 ±0.042) 


• 10^ 


(3.293 ± 0.043) 


• 10^ 


0.035 


(2.885 ± 0.037) 


• 10^ 


(2.553 ±0.039) 


• 10^ 


0.040 


(2.164 ±0.033) 


• 10^ 


(1.947 ± 0.034) 


• 10^ 


0.045 


(1.754 ±0.026) 


• 10^ 


(1.580 ±0.027) 


• 10^ 


0.050 


(1.314 ±0.025) 


• 10^ 


(1.202 ±0.025) 


• 10^ 


0.055 


(1.024 ±0.021) 


• 10^ 


(9.508 ±0.213) 


•10^ 


0.060 


(8.293 ±0.292) 


•10^ 


(7.692 ± 0.296) 


•10^ 


0.065 


(6.307 ± 0.300) 


■lOi 


(5.945 ± 0.304) 


•10^ 


0.070 


(4.636 ±0.180) 


■lOi 


(4.445 ±0.184) 


•10^ 


0.075 


(3.516 ±0.117) 


■10^ 


(3.430 ±0.114) 


•10^ 


0.080 


(2.673 ±0.115) 


■10^ 


(2.560 ±0.110) 


•10^ 


0.085 


(2.271 ±0.216) 


■10^ 


(2.214 ±0.217) 


•10^ 


0.090 


(1.412 ±0.204) 


•10^ 


(1.395 ±0.206) 


•10^ 


0.095 


(1.085 ±0.057) 


•lOi 


(1.056 ±0.059) 


•10^ 


0.100 


(7.412 ±0.584) 


■10° 


(7.377 ± 0.588) 


•10° 


0.105 


(5.069 ±0.537) 


•10° 


(5.121 ±0.529) 


•10° 


0.110 


(2.817 ±0.400) 


• 10° 


(2.914 ± 0.399) 


•10° 


0.115 


(2.652 ±0.329) 


•10° 


(2.428 ±0.301) 


•10° 


0.120 


(1.353 ±0.221) 


• 10° 


(1.516 ±0.183) 


• 10° 



We show the next-to-leading order perturbative predic- 
tion in QCD for S(j/4) in Fig. 7. In the same figure, the 
inset shows the K{y4) factors of the distributions. The 
physical cross sections for the two algorithms are very 
similar. The K{yi) factors are quite large, but much 
smaller than in the case of other four-jet event shape 
distributions. They depend weakly on the j/4 value for 
2/4 > 0.1 and decrease rapidly with decreasing 1/4 below 
2/4 = 0.1. In the case of the Cambridge algorithm the 
radiative corrections are 15-30% smaller than those for 
the Durham algorithm. 

In order to define Tminj we first have to define the 
thrust and thrust major axes The thrust axis ut is 
the direction ft which maximizes the expression 



T = max . ^1^1 



Eg \Pa 



(58) 
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FIG. 7. The next-to-leading order QCD prediction for the 
(solid) and j/f (dashed) differential distributions with 
renormalization scale x^j, — 1. The inset shows the K fac- 
tors of the distributions. 



where the sum runs over all final state hadrons (or par- 
tons). The thrust major axis is a three- vector nr^-^^j for 
which the expression in Eq. ( |5^ ) is maximal with the con- 
straint that nT,„aj is perpendicular to Ht, nr^^^j ■ Ht — 0- 
In order to obtain the value of T„iin, one evaluates the 
expression in the parentheses for a vector perpendicular 
to both fiTmaj ^^"i ^T- 

The C parameter is derived from the eigenvalues 
of the infrared safe momentum tensor 



PaPj 
\Pa\ I 



(59) 



where the sum on a runs over all final state hadrons and 
p\ is the ith component of the three-momentum pa of 
hadron a in the cm. system. The tensor is normalized 
to have unit trace. In terms of the eigenvalues \i of the 
3x3 matrix 9, the global shape parameter C is defined 
as 



C — 3 (A1A2 + A2A3 -f- A3A1) 



(60) 



The kinematical limit of the C parameter for three- 
parton processes is C — 0.75. Therefore, in the region 
C G [0.75, 1] the four-parton processes contribute to the 
leading order prediction, and our program is capable to 
calculate the radiative correction to the distribution. The 
results of such a calculation for the Born functions ST,„i„ 
and Bq agree with the known results (see e.g., |^^) The 
CT„,i„ and Cc correction functions are given in Table III. 



TABLE III. Correction functions to the differential distri- 
butions of the Tmin and C parameter event shape variables. 
The parameter values are at the lower edge of the correspond- 
ing histogram bin. 









C 


Cc 




0.00 






0.75 


(4.775 ± 1.100) 


■ 10^ 


0.02 


(3.319 ±0.270) 


■ 10'' 


0.76 


(6.082 ±0.160) 


■ 10^ 


0.04 


(2.381 ±0.082) 


■ 10'' 


0.77 


(4.610 ± 0.089) 


• 10^ 


0.06 


(1.652 ±0.038) 


■ 10" 


0.78 


(3.663 ± 0.063) 


• 10^ 


0.08 


(1.172 ±0.025) 


■ 10" 


0.79 


(2.904 ± 0.042) 


• 10^ 


0.10 


(8.600 ±0.130) 


■ 10^ 


0.80 


(2.406 ±0.031) 


• 10=^ 


0.12 


(6.488 ± 0.100) 


■ 10^ 


0.81 


(1.948 ±0.026) 


■ 10^ 


0.14 


(4.695 ± 0.077) 


• 10=* 


0.82 


(1.625 ±0.024) 


■ 10=* 


0.16 


(3.499 ± 0.042) 


• 10=* 


0.83 


(1.365 ±0.023) 


■ 10=* 


0.18 


(2.684 ± 0.027) 


■ 10^ 


0.84 


(1.135 ±0.019) 


• 10^ 


0.20 


(2.010 ±0.021) 


■ 10^ 


0.85 


(9.194 ±0.130) 


■ 10^ 


0.22 


(1.498 ±0.017) 


■ 10^ 


0.86 


(7.906 ±0.110) 


■ 10^ 


0.24 


(1.122 ±0.013) 


■ 10=^ 


0.87 


(6.293 ± 0.092) 


■ 10^ 


0.26 


(8.247 ± 0.100) 


■ 10^ 


0.88 


(5.217 ±0.084) 


■ 10^ 


0.28 


(6.093 ± 0.074) 


■ 10^ 


0.89 


(4.296 ± 0.066) 


■ 10^ 


0.30 


(4.501 ±0.180) 


■ 10^ 


0.90 


(3.391 ±0.052) 


• 10^ 


0.32 


(3.026 ± 0.057) 


■ 10^ 


0.91 


(2.815 ±0.061) 


• 10^ 


0.34 


(2.229 ±0.050) 


■ 10^ 


0.92 


(2.075 ± 0.057) 


■ 10^ 


0.36 


(1.549 ± 0.046) 


■ 10^ 


0.93 


(1.626 ± 0.032) 


■ 10^ 


0.38 


(1.095 ±0.028) 


■ 10^ 


0.94 


(1.221 ±0.026) 


■ 10^ 


0.40 


(7.100 ±0.210) 


■10^ 


0.95 


(8.154 ±0.260) 


■10^ 


0.42 


(4.437 ±0.180) 


■10^ 


0.96 


(5.193 ±0.190) 


•10^ 


0.44 


(2.684 ±0.190) 


■10^ 


0.97 


(3.165 ±0.130) 


•10^ 


0.46 


(1.439 ±0.150) 


•IQi 


0.98 


(1.312 ±0.094) 


■10^ 


0.48 


(6.447 ± 0.560) 


•10° 


0.99 


(2.769 ± 0.260) 


■10° 



In the case of event shape differential distributions the 
next-to-leading order corrections should logarithmically 
diverge at the edge of the phase space. This divergence 
occurs at zero for the 2/4 and Tmin distributions and is 
regularized by the multiplication with the value of the 
variable (see Eq. (|55|)). This is not the case for the C 
parameter, because it diverges at C = 0.75. Nevertheless, 
we obtained a finite and positive contribution in the first 
bin owing to bin smearing as we have checked explicitly 
by refining the bin width. 

Figs. 8 and 9 show the leading and next-to-leading or- 
der QCD prediction for the Tmin and C parameter dif- 
ferential distributions at = 1. The insets show the 
K factors which are large in both cases indicating 100 % 
or larger radiative corrections. As a result, the renor- 
malization scale dependence remains large, only the ab- 
solute normalization of the distributions increases with 
a factor of more than 2 with the inclusion of the radia- 
tive corrections. This feature is demonstrated in Fig. 10, 
where we show the scale dependence of the leading and 
next-to-leading order prediction for the average value of 
the thrust minor (above Tmin = 0.02) and C parameter 
(above C = 0.75). The leading and next-to-leading order 
curves run almost parallel down to ~ 0.2, only the 
latter is shifted to larger values. 
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FIG. 8. The leading-order (dashed) and the next-to-leading 
order (solid) QCD prediction for the Tmin variable with renor- 
malization scale a;^ = 1. The inset shows the K factor of the 
distribution. 
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FIG. 9. The leading-order (dashed) and the next-to-leading 
order (solid) QCD prediction for the C parameter with renor- 
malization scale = 1. The inset shows the K factor of the 
distribution. 
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FIG. 10. The renormalization scale dependence of the av- 
erage values of < Tmin >o.02 and < C >o.75 at leading and 
next-to-leading order. 



C. Radiative corrections to four-jet observables: 
summary 

In this subsection we summarize the resuhs of our ra- 
diative correction calculations for the various four-jet like 
distributions presented in previous publications and in 
this article. 

The QCD prediction at tree level (with renormaliza- 
tion scale = 1) is in general falls significantly be- 
low the measured values for unnormaliscd distributions 
of four-jet observables. Consequently, the calculation of 
the ncxt-to-leading order corrections to these cross sec- 
tions is indispensable for attempting a serious compari- 
son between data and theory. Our calculations show that 
the corrections are very large and the agreement in the 
comparison improves considerably with the inclusion of 
the radiative corrections. In particular, we found the fol- 
lowing features: 

(1) In the case of four-jet rates, the radiative corrections 
are about 100% for JADE- type clustering algorithms 
]ll| , p^ , while for the Durham algorithm it is less than 
60 % and even smaller for the Cambridge algorithm. The 
scale dependence for the latter algorithms is substantially 
reduced. The agreement between data and theory for the 
Durham clustering is very good and extends to small val- 
ues of j/cut when one matches the fixed order prediction 
with improved resummed next-to-leading logarithmic ap- 
proximation. 

(2) In the case of event shape variables the corrections 
are usually more than 100 % (the K factors are larger 
than 2). The residual renormalization scale dependence 
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is large indicating that even higher orders are important. 
One may conclude that, with the exception of the jet- 
related 2/4 distributions, these distributions cannot be re- 
liably calculated in fixed order perturbation theory and 
cannot be used for precision tests of QCD. 
(3) In the case of normalized angular distributions the 
corrections are small as expected (the K factors are close 
to 1) p^ , p^ . The renormahzation scale dependence is 
small, which however, does not mean that the effect of 
the radiative corrections on the measurement of the QCD 
color charges is negligible. According to Ref. |l^, the 
measured value of the Tr/Cf ratio may differ up to 25 % 
when leading, or next-to-leading order QCD predictions 
are used in the color charge fits. 

IV. CONCLUSIONS 

This paper dealt with the next-to-leading order cal- 
culation of four-jet observables in electron-positron an- 
nihilation. We gave details of the analytical calculation 
that lead to the construction of a Monte Carlo program 
PI which can be used to calculate the differential distri- 
bution of any four-jet observable at the 0(q;^) accuracy. 
The dipole method was used for achieving the analytical 
cancellation of infrared divergences. We described that 
modification of the algorithm which we found useful from 
numerical point of view. However, the modification is not 
essential as far as the theory is concerned. 

Compact formulas were presented for the Born-level 
five-parton helicity amplitudes and for the Born-level 
four-parton color-correlated matrix elements which are 
necessary for other next-to-leading order calculations, 
such as the next-to-leading order cross section of three-jet 
production in deep inelastic scattering and that of vec- 
tor boson plus two-jet production in hadron collisions. 
We also gave a group independent decomposition of the 
Born-level five-parton matrix elements. 

We calculated the next-to-leading order corrections to 
the four-jet rates with the Durham and Cambridge jet 
clustering algorithms and to the differential distributions 
of the 2/4, thrust minor and C parameter (at C > 0.75) 
four-jet shape variables. In the case of four-jet rates the 
radiative corrections were found to be large, but just 
acceptable. The renormahzation scale dependence de- 
creased significantly and the fixed order result matched 
with the next-to-leading logarithmic approximation gave 
remarkably good agreement with LEP data over a wide 
range of the resolution variable. The high level of agree- 
ment implies that the QCD four-jet background to 
pair production at higher center of mass energies can be 
predicted in perturbation theory reliably. In the case of 
event shape variables the radiative corrections and the 
renormahzation scale dependence are unacceptably large 
suggesting that the next-to-leading order prediction is 
not reliable and even higher orders are important. 
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APPENDIX A: HELICITY AMPLITUDES 

In this appendix we present analytic formulas for the 
four- and five-parton tree-level helicity amplitudes of the 
relevant subprocesses. These amplitudes were first cal- 
culated in Refs. The reason for presenting our re- 
sults here is twofold. On one hand we express the rele- 
vant color subamplitudes in terms of Weyl spinors |A;±), 
which were also employed in the case of the one- loop four- 
parton amplitudes ]20| , while on the other we found that 
our expressions in the case of the four-quark processes 
are more compact and the corresponding computer code 
is faster than earlier ones. Another new feature of the 
amplitudes in this appendix is that we allow for the ex- 
istence of light fermionic degrees of freedom in the ad- 
joint representation of the gauge group (light gluinos). 
In calculating the amplitudes, we used quark and gluon 
currents and standard helicity techniques p3| , p^ . 

We consider three subprocesses, each involving a vec- 
tor boson V{Q) carrying total four-momentum Q and n 
QCD partons (n = 4, or 5 here). The first subprocess 
is the production of a quark-antiquark pair and n — 2 
gluons. The second one is the production of two quark- 
antiquark pairs (of equal, or unequal flavor) and n — 4 
gluons. Finally, the third process is the production of a 
quark-antiquark pair, a light-gluino pair and n— 4 gluons: 

(+{-Pi) + t-{-Pl)^V{Q)^ (Al) 

qiVl) + .91 (P2) + . . • + gn-2{Pn-l) + q{Pn) , 

l+{-pe) + r{-pi) V{Q) ^ (A2) 
q{pi) + q{p2) + Q{p:i) + Q{Pi) 

+9l{Pb) + • ■ . +5„-4(Pn) , 

i+{-Pi)+r{-pj)^V(Q)^ (A3) 

q{pi) + q{p2) + gips) + givd 

+9l{Pb) + • ■ ■+9n-i{Pn) ■ 

We have chosen the crossing invariant all particle outgo- 
ing kinematics with corresponding particle-antiparticle 
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assignment, therefore, momentum conservation means 

Pe+Pe+Pl+P2+P3+P4:+P5 + ■■■+Pn = ■ (A4) 

We shall express the amplitudes in terms of color sub- 
amplitudes. In the case of process (A_l), the color basis 
is chosen to be product of generators in the fundamental 
representation (in this appendix we use the normaliza- 
tion Tr = 1 in Tr(rt^) = TrS'''' for the generators of 
the symmetry group), therefore, the helicity amplitudes 
have the decomposition: 



Na 



I -Lj 1 J • ■ • J 



(A5) 



^ (t-.....r"-),^,„m(l^S...,nJ"), 

{2,...,Tl-l} 

where {2, . . . , ti — 1} denotes all permutations of the la- 
bels (2, . . . , n — 1) and m(l, . . . , n) are the color subam- 
plitudes. In Eq. (A5) and in the following formulas the 
lepton labels are suppressed. 

In the case of the four-fermion subprocesses (processes 
( [A^ ) and (A3)) we decompose the helicity amplitudes as 
follows: 



(A6) 

E E(-l)'' E(-l)''-^"(l'2'3,4,5,...,n). 

{5,....ri} {1,3} {2,4} 

where P = if the elements are in the canonical order 
((1,3), or (2,4)) and P — 1 if the elements are permuted 
((3,1), or (4,2)). The partial amplitudes An can be de- 
composed further in color space. In the case of four-quark 
production, 

A(1„2,-,3q,4q) = (A7) 

where M(lJ^\ 2^^" , 3^^^ 4^^*) are the color subamplitudes 
and T(l, 2, 3^ 4) fs defined by 



Na 



T(l,2,3,4) = ^<,^i 



b 

13 U 



(A8) 



6=1 



In the case of four-quark plus one-gluon production, there 
are four independent basis vectors in color space: 



A(l<?,2g,3Q,4Q,5g) = 

4 

r,(l, 2, 3, 4, 5) M,{\)l , 2)1 , 3^3^ , 4^^, 5''^) 



(A9) 



where Ti(l,2,3,4,5) are the color basis vectors: 



Na 



Ti(l,2,3,4,5) = ^(t''^i''),,,j, t: 



b _ 



(AlO) 



6=1 



6 _ 
1314 ' 



T2(l,2,3,4,5) = ^(t''i'^^),^^, t: 
6=1 

Na 

T3(l,2,3,4,5) = ^<,^ {eH\^j^, 
6=1 

Na 

T4(l,2,3,4,5) = 5]<,^ {th'^^U^. 



(All) 
(A12) 
(A13) 



6=1 



The partial amplitudes for the process (A3) can be 
written in terms of the color subamplitudes of the process 
(|A2|), only the color basis differs. When n = 4, 



-^4(lq,2g, 3g,4g) 



(A14) 



T(l, 2, 3, 4) M{l)l , 2)1 , 3''=' , 4''* 



where 



Na 



T(1,2,3,4) = ^<,^P,\ 



(A15) 



6=1 



Finally, for n — b we have 

A(l?,2,-,3g,4g,5g) = 

4 



(A16) 



T^{l, 2, 3, 4, 5) M,{l)l , 2^;^% 3''% 4'^% 5"=) , 



where 



Na 



n{l,2,3,4,5)^Y(*-"^%hF. 
6=1 

Na 

r2(l,2,3,4,5) = 5](iV^),^j,F, 



6 

0304 ' 



6 

0304 ' 



6=1 
Na 



r3(l,2,3,4,5) = ^<.^ (F'^'^F'), 
6=1 

Na 

Ul,2,3A,5) = J2tU, {F'F'^-), 



(A17) 
(A18) 
(A19) 
(A20) 



6=1 



In the following subsections we give explicit formulas 
for the color subamplitudes with a common coefficient 
factored out: 



(A21) 



Mil)l,2J^''\3'}l4j!'')= (A22) 

M,{l)l , 2-''i , , 4-''^' (A23) 
2e^g^cYJl' %/4 -^^(l"^2""^3''^4""^5''^) , 



13 



(A24) 



M,(lJ^^27^''^3"^4-"^5"=) = (A25) 

2e2g3c}';)^'^ - A,(l''\2~''i,3''%4-''%5''^) , 



with s ^ ^ {pe +pi)'^. The C coefficients contain the 
electroweak couphngs. If the vector boson y is 7 or 
this coefficient is defined by 



= {-Q^''+v1'v';{'Vz{s))Sf,f,, (A26) 



where /i , /2 are the flavour indices of the quark antiquark 
pair that couples to the vector boson and 



-l + 2sin2 e 



w 



2 sin^ 9w 
sin 29w 

2Qf sin^ e 



(A27) 



w 



' sui2ew ' ' 

±1 - 2gf sin^eiw , 

Vf = , Vf — 

^ sin26'vF sin26w 

are the left- and right-handed couplings of leptons and 

3y/ 



quarks to neutral gauge bosons. In Eqs. (A27 A28) 9 



denotes the Weinberg angle, Q/ is the electric charge 
of the quark of flavor / in units of e and the two signs 



in Eq. (A28) correspond to up and down (— ) type 



quarks. The coupling C contains the ratio of the and 
photon propagators, 



s-Mj+ iTz Mz ' 



(A29) 



where Mz and Tz are the mass and width of the . 

If the vector boson is a or a W~ , then the 
couplings take the form 



(A30) 



1. Four-parton color subamplitudes 



In this subsection, we present all four-parton color sub- 
amplitudes for the helicity configuration hq — + and 
hi = +. The amplitudes for the reversed helicity config- 
urations can be obtained from these amplitudes by ap- 
plying parity operation P, which amounts to making the 
substitutions {ij) = {k~\kj') ^ [ji] = {kj'\k~). The am- 
plitudes when only the lepton helicities are reversed can 
be obtained simply by exchanging the lepton labels and 
flipping the lepton helicity in the coupling factors C'/j*}^^^ 
We use the notation 

{z\lm . . . \j) ^ C(fcn7M7. • ■ ■ 1^,^) , (A33) 



[t\lm...\j]^k^k-;^{k+\j^j,...\kf), 



(A34) 



m + m) . . . \j) EE (A:f + fc^^J . . . (fcr I 7m ■ • • \kf) , (A35) 

m + m) . . . \j] EE (fcf + k^J . . . (fc+l 7m . . . \kf) , (A36) 
(A28) and the two- and three-particle invariants s.y = {ki + kjy 



and tiji = {ki + kj + kiY . Labels 5 and 6 refer to the 
positron and electron respectively. 

The two-quark two-gluon amplitudes are as follows: 

(45) 2 [56] 



^(1^2+ 3+,4^) = - 



(12)(23)(34) 



2g , 3g , 4^ ) 



[16]2(56) 
[12] [23] [34] 



^(1(J^' 2g , 3g , 4^ ) — — 



(31)[12](45)(3|(1 + 2)|6) 



(A37) 



(A38) 



(A39) 



(12)s23/:i23 

(34)[42][16](5|(3 + 4)|2) , (5|(3 + 4)|2)(3|(1 + 2)|6) 



[34] 523^234 



(12)[34]S23 



where /i denotes the partner of quark /i in the 
SU(2)i doublet and, for the sake of simplicity, we set 
the Koba yashi-Maskawa mixing matrix to unity. In 
Eq. ( A30 ) the left- and right-handed co uplings di ffer from 
the corresponding expressions in Eqs. (A27,A2£): 



2V2sin6' 



w 



w+ = w+ = . 



(A31) 



In this case Vw{s) denotes the ratio of the and pho- 
ton propagators. 



Vwi-s) = 



s - + iTw Mw 



(A32) 



where Mw and Fw are the mass and width of the W^. 



^(Ig^' 2g , 3g , 4^ 



[13]2(45>(21(1 + 3)16) 



[12]s23i 



123 



(A40) 



(24)2[16](51(2 + 4)13) [13] (24) [16] (45) 

(34)523*234 [12](34)S23 

The four-quark amplitudes are as follows: 

A(l+,2,,3+,4^)= (A41) 

Al3](52)(4l(l + 3)16) , (42)[61](51(2 + 4)13)\ 



\ il34S34 



^234534 



2g , 3q, 4^) — A(l+, 2^ , 4(^, 3q) . 



(A42) 
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2. Five-parton color subamplitudes 



In this subsection, we present all five-parton color sub- 
amplitudes for the helicity configuration hq = + and 
hi = +. The amplitudes for the remaining helicity con- 
figurations can be obtained from these amplitudes as in 
the n = 4 case. Labels 6 and 7 refer to the positron and 
electron respectively. 

First we list the two-quark three- gluon amplitudes: 
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2+ 3+ 4+ 5,-)- ^12) (23) (34) (45) ' ^"^^^^ 

^(1^2^3^4;,5r) = tyj(f+,?l^) (mm±m + Mi+M) (A44) 
+ 



(23)(34)[34]t567 V i234 (12) 

(4|(5 + 6)|7)(6|(5 + 4)|3) (4|(1 + 2)|7)(6|(5 + 4)|3)(45)[53] [17](6|(1 + 7)|2)(45)^[53F 
(12) (23) (34) [34] [45] (12)(24)[45]s34i345 [45] (42)534*345^167 



, [1"](^5)[.53] / (Gl(l + 7)|2) (6|(1 + 7)|3) \ _ [17] (45) [23] / (G|(l + 7)|2)(24) 
^(23)[34][45]ti67 I (34) ^ (24) ^ (23)[34]i234ii67 I (34) 



+ + + . ^ (31)[12](3|(1 + 2)|4)(3|(5 + 6)|7)(65) (31)[12](3|(1 + 2)|7)(65)(35) 

^' ^' ^' (12)(34)S23W567 (12) (34) (45)523^123 



(3|(1 + 2)|7)(6|(1 + 7)|2)(35)^ (3|(1 + 2)|7)(6|(5 + 3)|4)(35)[42] _ [42]^[17](6|(1 + 7)(2 + 4)|3)(35) 

(12)(23)(34)(45)[23]t345 [23](12)(23)s34i345 523534^234^167 

^ [42](3|(5 + 6)|7)(65) / [42](3|(2 + 4)|1) _ (3|(1 + 2)|4) \ ^ [17](6|(1 + 7)|2)(35)^ / (3|(5 + 4)|2) ^ [42][54] \ 

S23S34t567 V ^234 (12) / 523*345*167 V (34) (45) S34 / 



,.^+ „_ ,+ _ [43]^(2|(3 + 4)|1)(2|(5 + 6)|7)(65) [13](2|(3 + 4)|1)(2|(5 + 6)|7)(65) 

^(l,'^S"^.'4,,5,j- 523534*234*567 ^ [12] (34) (42)^23*567 ^ ^ 



[13]^(2|(1 + 3j|4)(2|(5 + (.i)|7)((i5) [13]^(2|(1 + 3) 1 7) (65) (25) [17] (25) [43]^((.| (1 + 7)i3 + 4j|2) 

[12] (24)523*123*567 [12] (24) (45)523*123 523^34*234*167 



(61(5 + 4)13)^-11 



(42) 

[12]523(24)t345 V'""" ' V(45) (34); (6|(5 + 3)|4)-py 

^[17](6|(1 + 7)|3)(25) A ^ ^^|3^ f " gl V (2|(5 + 3)|4) ^^^^ 



(24)523*345*167 V V (45) (34) ^ ^'^ ^''(34) 

2+ 3- 4- 5-^ = [12][2|(3 + 4)(5 + 6)|7](65) / (43)^ ^ (31) \ _ (31)^[12]^(4|(5 + 6)|7)(65) 

^ ^' 523*567 ^534*234 [34] [42] (12) ^ (12) [24]523*123*567 ^ ' 

(31)[12](3|(1 + 2)|7)(6|(5 + 4)|2) (3|(1 + 2)|7)(6|(1 + 7)|2) [17](6|(1 + 7)|2)(3|(5 + 4)|2) 
(12)[24][45]S23*123 (12)[34][45]S23 [34] [45] 523*167 

[17](6|(1 + 7)|2)(5|(3 + 4)|2)(43)^ ^ 

523534*234*167 

523534*234*567 

[14](3|(5 + 6)|7)(65) /[42](2|(1 + 3)|4) + [43](3|(1 + 2)|4) [14](3|(1 + 2)|4) 



[42]S34*123*567 V [23] [12] 

[14](65)(35) /[42](2|(l + 3)|7) + [43](3l(l + 2)]7) [14](3|(1 + 2)|7) \ , [14][17](6|(5 + 3)|4)(35)2 



[42] (45)534*123 V [23] [12] / [12] [24] (45)534*345 

[17](6[(l + 7)[4)(2|(5 + 3)[4)(35)^ _ [17](6[(1 + 7)[4)(5[(2 + 3)[4)(35) _ [17](6[(1 + 7)[4)(5[(2 + 3)[4)(23)^ 

[42] (45)534*345*167 (45) [42] [23]534*167 523534*234*167 
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^(1+ 2;,3+ 4;,5,) = PA(5+,4+ 3^,2+ ir) , 2^, 3^-, 4^ , 5, ) = P A(5+, 4+, 3+ 2+ ir) . (A49) 

|6<->7 |6<^7 

The four-quark one-gluon amplitudes have the form: 

, „_ „+ ,_ _ [15](4|(1 + 5)|3)(4|(2 + 6)|7)(62) [17](6|(1 + 7)|5)(42)^[23] 

(4|(1 + 5)|7)(6|(2 + 4)|3)(42) ^ [53](4|(3 + 5)|1)(4|(2 + 6)|7)(62) ^ [17](6|(1 + 7)(3 + 5)|4)[35](42) ^ 

(15)(54)s34t234 (45)334^345^267 (45)s34t345il67 

„+ ,_ _ [13]^(51)(4|(2 + 6)|7)(62) [17](6|(1 + 7)|3)(5|(2 + 4)|3)(42) 

' ' "^Q' 5. ) - [35]5i5S34i267 + [35]s34t234tl67 ^ ^ 

[13][17](6|(2 + 4)|3)(42) ^ [13](54)[3|(4 + 5)(2 + 6)|7](62) [17](6|(1 + 7)|3)(54)(2|(4 + 5)|3) ^ 

[15][53]s34t234 [35] 534*345^267 [35] 534*345*167 

2r, 3q, 4+5+) = Ai(l+ 2r, 4+ 3^, 5+) , ^i(l+ 2,-, 3q, 4+, 5;) = Ai(l+ 2^, 4+ 3^, 5;) , (A52) 
, ,_ _ [53](4|(3 + 5)|1)(4|(2 + 6)|7)(62) [17](6|(1 + 7)|3)[25](42)^ 

^2(1,,2,,3q,4^,5J- (45)534*345*267 (45)525^34*167 ^ ' 

[13](4|(1 + 3)|5)(4|(2 + 6)|7)(62) [13](4|(1 + 3)|7)(42)(62) [17](6|(1 + 7)(3 + 5)|4)[35](42) ^ 

(45)534*134*267 (45) (52)534*134 (45)534*345*167 

[13](54)[3|(4 + 5)(2 + 6)|7](62) [13]^(41)(5|(2 + 6)|7)(62) 

^2(1,, 2,, 3q, 4^, 5,)- [35]534*345*267 ^ [35] 534*134*267 ^ ' 

^ [17](6|(1 + 7)|3)(54)(2|(4 + 5)|3) _ [13](4|(1 + 3)|7)(6|(2 + 5)|3) ^ [17](6|(1 + 7)|3)(52)(4|(2 + 5)|3) ^ 

[35] 534*345*167 [35] [52] 534*134 [35] 525 534*167 

^2(1+ 2r,3Q,4+ 5+)=^2(l+ 2r,4+ 3^,5+), ^2(1+ 2r, 3^, 4+, 5,") = ^2(1+ 2r, 4+ 3^, 5,") , (A55) 
A fl+ 2- 3+ 4- 5+) = [53](4|(3 + 5)|1)(4|(2 + 6)|7)(62) _ [17](6|(1 + 7)(3 + 5)|4)[35](42) 

'5' 9^ (45)535*345*267 (45)535*345*167 ' ^ ' 

, (,+ ^- 4+ .-X _ [14](53)[4|(3 + 5)(2 + 6)|7](62) [17](6|(1 + 7)|4)(53)(2|(3 + 5)|4) 

^3(l,,2,,3g,4^,5,)- [45]535*435*267 ^ [45]535*435*167 ' ^ ^ 

, „_ „+ ,_ _ [13] (54) [3| (4 + 5) (2 + 6) 1 7] (62) [17](6|(1 + 7)|3)(54)(2|(4 + 5)|3) 

^4(1,,2,,3q,4^,5J- [35]545*345*267 [35]545*345*167 ' ^ ^ 

, ,+ + [54](3|(4 + 5)|1)(3|(2 + 6)|7)(62) [17](6|(1 + 7)(4 + 5)|3)[45](32) 

^4(1, , 2, , 3q, 4^, 5 J - (35)545*345*267 ^ (35)545*345*167 ^ ^ 
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APPENDIX B: MATRIX ELEMENTS 



1. Four-parton color-summed matrix elements 



In this appendix we present analytic formulas for the 
color-correlated four-parton Born-level matrix elements 
and for the four-, five-parton Born-level matrix elements. 
The calculation of the color-correlated four-parton am- 
plitudes is a straightforward application of color alge- 
bra and the four-parton helicity amplitudes. However, 
to our knowledge these results were not published pre- 
viously. The uncorrelated color sum was first calculated 
in Ref . . We present our results in terms of the color 
subamplitudes of Appendix A. It is a new feature of the 
matrix elements in this appendix that they are given in 
terms of group independent functions and eigenvalues of 
the quadratic Casimir operators of the underlying gauge 
group. 

Having the helicity amplitudes at our disposal, we cal- 
culate the squared matrix elements summed over final 
state colors without and with color-correlation: 



\Mnil,...,n)\' 



((1, . . . . . . ,n)„ , n = 4,5. 



A<r(l,...,4) 



(Bl) 



(B2) 



1 'ii 



I T, ■ T, 



11"'^ 



where in the latter case we leave the helicity indices ex- 
plicit so that both correlated and uncorrelated helicity 
summation is possible. (Although we did not show the 
flavor indices, the flavor summation is also left out, as 
will become clear later.) In the correlated case we have 
to insert the helicity matrix (see Eq. (|^)) 



H 



hh' 



Jhih' 



{h^\V.^,\h';) 



and in the uncorrelated case 



H 



hh' 



^hih[ 



Jhih' 



Jhr,h' 



(B3) 



(B4) 



We evaluate the color sum in such a way that the ma- 
trix elements are given as polynomial expressions of the 
Casimir invariants of the gauge group with group inde- 
pendent kinematical coefficients. In addition to the usual 
quadratic Casimirs Cp and Ca, we shall also use a cubic 
Casimir C3 that is defined as 



Na 



Cs^ J2 Tr(t'^i''r)Tr(t'=tV) 



(B5) 



a,b,c—l 



In the following subsections we list the explicit formulas 
for \M4\\ \M'i\^ and {M^l^. 



In this subsection, we give explicit formulas for the 
color-summed Born matrix elements for four final state 
partons. There are four different cases: the two-quark 
two-giuon process and three four-fermion processes (two 
unequal flavor quark pairs, two equal flavor quark pairs 
and the two-quark two-gluino production). The color 
summation is straightforward in each cases, we simply 
list the results: 

\Miilq,2g,3gA<if ^ N,Cl (B6) 
x{|m(l;,,2,3,4/J|2 + |m(ly,,3,2,4/J|2 
+2 Re (m(l/, , 2, 3, 4/Jm(l/, , 3, 2, 4/J*) 
~x Re(m(l/, , 2, 3, 4/Jm(l/, , 3, 2, 4/J*) } , 



|A^4(lg,2,-,3Q,4Q)|'=iVcC| (B7) 
X { - 2 Re(M(l/, , 2/„ 3/3, 4/jM(l/, , 4^,, 3^3 , 2;,)*) 
+x Re(M(l/, , 2f, , 3/3 , 4/jM(l/, , 4/, , 3/, , 2/J*) 
+y (|M(l/,,2/,,3/3,4/J|' 

+ |M(l/,,4/„3/3,2/j|')}, 



|A^4(lg, 2,-, 3^,43)1 



(B8) 



M(l/,,2/,,3,4) 



where x and y are ratios of the quadratic Casimirs (see 
Eq. H) and M(l/i , 2/2 , 3/3 , 4/J is defined by 

M(l/,,2/„3/3,4/J= (B9) 
, 2/, , 3/3 , 4/J + M(3/3 , 4/„ 1/, , 2/, ) . 



2. Four-parton color-correlated, color-summed 
matrix elements 



In this subsection, we give explicit formulas for the 
color-correlated, color-summed Born matrix elements for 
four final state partons. We consider those four cases 
as in the previous subsection. The color summation is 
again fairly straightforward, therefore, we only record the 
results. 

For the V qqgg subprocess 



Xr(l,,2„33,4,)| = 



(BIO) 
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where the non- vanishing elements of the matrices Mq' , 

MJ^, M^^ arc given by 





= 2(5i +S2+ 53) , 


(Bll) 


7\ /rl2 
M^ 

Ml' 


— — ^Oi — ZO2 — 0O3 , 

= M^" = Mf = Mf = 5i + 52 + 53 


(B12) 


M^^ 


= -{S1 + S2 + 2S3) , 


(B13) 


XX 


= Mj^ = -i(52 + 53), 




XX 


= m!1 = ~{s, + Ss), 




XX 


= l{Si + S2) 





and the Si functions are defined by 

5i =m(1^^^2''^3''^4^^^)*m(l;;,2'^^3'^^4;;) , (B14) 

= m(1^^^3''^2''^4J^^rm(l;;,3''^2'^^4;^) , (B15) 

S3 = m{l'}l,2'^\3''^A'}:rm{l'}l,3'''^XK4:) (B16) 
+m{l'}l , 3''^ , 2^2 , 4^^^ )*m(l^^ , 2^= , 3^*3 , 4^^^ ) . 

In the case of four-quark production 

|Mf(l„2„3Q,4Q)|'= (B17) 



ik 



WMg+xyMll}, 



where the non-zero element of the matrices Mq'', MJ^ , 
Mi'', Mf , Mil, Mil following ones: 



M^' 
Ml' 




= ^ 

= 2^3, 


M',' = 


-2^3, 




(B18) 


Ml' 
Ml' 


= Mf 


= Ml' = 
= -3^3 , 


Mf = 


2^3, 




(B19) 


y 


= Mf 


= 25i, 




= Mf 


= 252, 


(B20) 


Ml' 
Ml' 


= Mf 
= Mf 


= 2S'2 , 

= -2(5i + ^2) , 


= Mf 


= 25i, 


(B21) 


XX 


= M" 

XX 


= M23 = 

XX 

= ^3, 


= 

XX 


-\S3, 




(B22) 



= Mfy = -S, , Mil = ^11 = , (B23) 
Mll=Mll = Si + S2. 

For this case the 5, functions are defined as follows: 



S, = M{l)l , , 3^ , 4^ )*M(l5; , 2;^ 3?:, 4?; ) , (B24) 

52 = M{l)l,4!}l,?,)l,2)iyM{lf^,A%2,'}l,2f^) , (B25) 

53 = M{\)1 , 2J^^ 3^^, 4^;)*M(lj; , 4;j, 3;j, 2j^^) (B26) 

+M{i'fi , 4): , , 2)^= )*M(i;; , 2;^ , 353 ,4fj, 



where 



M(l/,,2/„3/3,4/J= (B27) 
M{lf, , 2/, , 3/3 , 4/J + M{lf, , 4/„ 3/3 , 2/J . 



Finally, for the V — > qqgg subprocess 



^3 



(B28) 



where the non-vanishing elements of the matrices M^'^ 
and M^^ are given by 



Mi^ = 2 



M(l/„2/„3/3,4/J 



(B29) 



M^^ = -M-"" = - 

XX XX 



M(l/,,2/„3/3,4/,; 
1 



M^^ = M^l = Mil = Mil = 



(B30) 
2 



M(l/„2/„3/3,4/J 



3. Five-parton color-summed matrix elements 

In this subsection, we give explicit formulas for the 
color-summed Born matrix elements for five final state 
partons. There are again four different cases: the two- 
quark three-gluon process, the production of two equal, 
or unequal flavor quark pairs plus a gluon and the two- 
quark two-gluino one-gluon production. 

In the case of the two-quark three-gluon process the 
color summation is straightforward and leads to the fol- 
lowing expression: 



|A45(l„2„3„4g,5,-)r = 



(B31) 



N^Cl \Mo--{Mi + 2Mo) + —{Mo + Mi + M2) ^ , 



where 
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Mo 



^ m(l/„2,3,4,5/,) 

{2,3,4} 



M2= |m(l/,,2,3,4,5/,)r 

{2,3,4} 



(B32) 
(B33) 



and 
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Ml = -2 M2 



(B34) 



-2 Re ^ Jm(l/,,2,3,4,5/,)*[m(l/,,2,4,3,5/J + m(l/,,3,2,4,5/J-m(l/,,4,3,2,5/j') I 

{2,3,4}' I ^ ^ J 



with {2,3,4}' denoting the cychc permutations of the 
three labels 2, 3 and 4. 

In the case of the four-quark onc-gluon subprocesses 
we have to evaluate the following color sums: 



tIti{2 ^ 4) = rjr2(i ^ 3) = 



N.C'p (Cf--Ca] , 

t/Ti(1 ^ 3) = T2^T2(2 ^ 4) = 

N,Cf (cf - ^Ca) {Cf - Ca) , 



(B35) 

(B36) 

(B37) 



t/Ti(1 ^ 3, 2 ^ 4) = T^T2{1 ^ 3, 2 ^ 4) = (B38) 



C3 - N.CfTr^ , 



tIt.2 = N.CfTr I^Cf - -CaJ , 

tIT2{1 ^ 3) = tIT2{2 ^ 4) = 
1 



N,Cf [ Cf - -Ca 



T|r2(1^3,2^4) =C3 



(B39) 
(B40) 

(B41) 



Using these results, the square of the matrix element for 
any flavor configuration can be written in the form: 



|7W5(lg,2,-,3Q,4Q,5g^ 



(B42) 

N^C% {Mo + xM^ + yMy + zM^ + x^M^^ + xyM^y} , 

where we have introduced the ratio 

C3 



N,Cl 

and the following abbreviations: 
Mo = B + C + E , 



M^ = -]^{'iC + 2E + B) 



My=A + D 



M, = F + G., 



(B43) 

(B44) 
(B45) 
(B46) 
(B47) 



M,, = -{2C + E), 



-AF + D), 



with the functions A, B, C, D, E, F defined as 

{1,3} {2,4} i=l 



(B48) 
(B49) 

(B50) 



B = -2Re(MiMi(2 ^ 4)* + M2M2(1 ^ 3)* (B51) 
+(1^3,2^4)) , 

C = -2Re(MiMi(l ^ 3)* + M2M2(2 ^ 4)* (B52) 
+(1^3,2^4)) , 



D = 2Re(^Y^ M1M2) , 

{1.3} {2,4} 

E = -2Re((Mi + Mi(l ^ 3, 2 ^ 4)) 
x(M2(l ^ 3) + M2(2 ^ 4))* 
+(Mi ^ M2)) , 

F = 2Re(MiMi(l ^ 3,2 ^ 4)* 
+Mi(l ^ 3)Mi(2 ^ 4)* 
+(Mi ^ M2)) , 

G = 2Re(MiM2(l ^3,2^ 4)* 

+Mi(l ^ 3)M2(2 ^ 4)* 
+(Mi ^ M2)) , 

where 



(B53) 
(B54) 

(B55) 

(B56) 



(B57) 



^i(l/i'2/„3/3,4j„5g) = 
_ M,{lf„2f„ 3/3 Au,5g) + M3 (3/3 , 4/„ 1/, , 2/, ,5g) , 
^2(1/1, 2/,, 3/3, 4/,, Sg) = 

^2(1/1, 2/, , 3/3 , 4/„ 5g) + M4 (3/3 , 4/, , 1/, , 2/, , 5g) . 

For the V qqfjfl!) process wc have to calculate the 
following products of the Ti colour factors: 

(B58) 
(B59) 





= tIt2 


= N.CICa 






= N.CfCI 
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TIT2 = N.CfCa { Cf - ^Ca 



TlTi = tITz = N^Cf 



tIt^=N,Cf'^. 



(B60) 
(B61) 
(B62) 
(B63) 



Using these identities the square of the matrix element 
can be written in the form: 



|Al5(l„2,-,3§,4g,5,)r 



(B64) 



where 



Ml + M2 



M3 + M4 

;Re((Mi + M2)(M3 + M4)* 



(B65) 
(B66) 



-Re((Mi+M4)(M2+M3)*) . 
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